function[Rhxy,Rhyx,Phxy,Phyx] = getrho1frq(tf,frq),

for i = 1:length(tf),   
   Phxy(i) = atan2(imag(tf(i,1)),real(tf(i,1))) * 180 / 3.1415926535897;
   Rhxy(i) = 0.2 * (1/frq) * (tf(i,1) * conj(tf(i,1)));
   Phyx(i) = (atan2(imag(tf(i,2)),real(tf(i,2))) * 180 / 3.1415926535897);
   Rhyx(i) = 0.2 * (1/frq) * (tf(i,2) * conj(tf(i,2)));
end;